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We propose a new mechanism leading to scale-free networks which is based on the pres- 
ence of an intrinsic character of a vertex called fitness. In our model, a vertex i is assigned 
a fitness Xj, drawn from a given probability distribution function f(x). During network 
evolution, with rate p we add a vertex j of fitness Xj and connect to an existing vertex 
i of fitness X{ selected preferentially to a linking probability function g(xi,xj) which de- 
pends on the fitnesses of the two vertices involved and, with rate 1 — p we create an edge 
between two already existed vertices with fitnesses Xj and Xj, with a probability also 
preferential to the connection function g(xi, xj). For the proper choice of g, the resulting 
networks have generalized power laws, irrespective of the fitness distribution of vertices. 
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Complex networks are powerful tools to describe a large variety of biological, 
social, and technical networks. A network is a mathematical object which consists 
of vertices connected by edges. Despite differences in their nature, many real- world 
networks are characterized by similar topological properties, in contrast to those 
obtained by traditional random graphs. One of the most interesting phenomena is 
the scale-free (SF) behavior, which means a power-law distribution of connectivity, 
P(k) ~ fe~T J where P(k) is the probability that a vertex in the network is of degree 
k and 7 is a positive real number determined by the given network. In order to 
understand how SF networks arise, much work has been done in the past decade. 
It has been shown that growth and preference seem to be the principal mechanisms 
for SF behavior. 

The exploring the preference can be directed in two classes. The first 
class of research is based on the rich- get-richer rule, which was implemented 
by newcomers preferent ial connecting to old vertices with certain topological 
characteristics! 1 l 2 l a l 4 l 5 El i n the best known Barabasi- Albert (B A) ModelP the net- 
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work grows at a constant rate and new vertices attach to old ones with probability 
H(i) ~ fcj. In this way, vertices of high degree are more likely to receive further 
edges from newcomers. In fact, this extreme assumption is not always available for 
many networks when their sizes are huge. The similar mechanism was also used 
in weighted networks, for instance, the network grows at a constant rate and new 
vertices attach to old ones with probability II(z) ~ s,, where Sj = ^2jeU(i) Wi i 
is the strength of vertex i and the sum runs over the set U(i) of neighbors of 

The second class of research utilizes the fit-get-richer mechanism, which was 
carried out by newcom ers preferential connecting to old vertices with high intrinsic 
fitnesses F^9 | 10 | ll | 12 |'j'his is better adapted to model certain networks where topo- 
logical properties are essentially determined by "physical" information intrinsically 
related to the role played by each vertex in the network, such as the ability of an 
individual, the content of a web page, or the innovation of a scientific article. 

Caldarelli et al recently introduced a varying vertex fitness modelj^ where they 
consider an undirected graph of N vertices. At every vertex i a fitness Xi, which is 
a real number measuring its importance or rank, is assigned. Fitnesses are random 
numbers taken from a given probability distribution f(x). For every couple of ver- 
tices, i and j, an edge is created with probability g(xi,Xj) (a symmetric function 
of its arguments) depending on the "importance" of both vertices, i.e., on a:,; and 
Xj . Actually this is a natural generalization of the classic Erdds-Renyi graphic Al- 
though it is a static model, the network recovers the power-law behavior of degree, 
betweenness, and clustering coefncientP On the other hand, Bedogne and Rodgers 
proposed a growing network with intrinsic vertex fitnesses!^ Besides employing 
the edge-created mechanism suggested in Ref. they also considered two cases of 
new vertices connecting to old ones, uniform or degree-preferential. The interplay 
between the fitness linking mechanism and uniform attachment results in an ex- 
ponential degree distribution for any fixed fitness x, while the degree- prefere ntia l 
attachment instead induces that the degree distribution decays as a power lawP^ 

Models of the first class often present us such an evolution picture: old vertices 
are passively attached by newcomers according to the degree- (strength-) preferen- 
tial mechanism. On the contrary, models belongs to the second class pay much at- 
tention to the creation and reinforcement of internal connections. Combining above 
two aspects, we argue that the connection between two vertices is the result of their 
mutual affinity and attachment. Not only for interactions among new vertices and 
old ones, but also for that among old vertices, which we call "mutual selection" . Mo- 
tivated by this, we suggest an evolving network model ruled by the fitness-dependent 
selection dynamics. The generated network has a good right-skewed distribution of 
degrees. 

The present model starts from an initial m isolated seeds and each vertex i is 
endowed with a fitness Xi > 0, drawn from a given probability distribution f(x). 
At each time step, we perform either of the following two operations, (i) With 
rate p € (0,1) we add a new vertex j of fitness Xj £ f(x) to the network. The 
new vertex connects to an existing vertex i of fitness x% selected preferentially to 
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a linking probability function g(xi,Xj) which is symmetric and dependent on the 
associativity of the both vertices, (ii) With rate 1 — p we create an edge between 
two vertices, i and j, already presented in the network with the probability also 
preferential to their integration g(xi,Xj). After t time steps, this scheme generates 
a network of m + pt vertices and t links. Notice that either process is chosen in the 
network growth, only one edge is added to the system at each time step (duplicate 
and self-connected edges are forbidden), however, this is not essential. 

In our model, each vertex is assigned a fitness, either initial seeds or subsequent 
newcomers. Denoting N k (x,t) the average number of vertices with degree k and 
fitness x at time t, we can write out the rate equation for network evolution 

dN k (x,t) = P J °° f(x')g(x,x')[N k ^{x,t) - N k (x,t)]dx' 
9t J °° f{x') £r=o Jo°° g(x, x')N k (x, t)dxdx> 

+ pSk,if(x) 

+ 2(l-p)[N k -x(x,t)-N k (x,t)] 



Y,T= Jo°° 9( x > x')N k {x', t)dx' 



Efclo Jo Efclo Jo°° 9 (x, x')N k {x', t)dx'N k (x, t)dx ' 



(1) 



The first term on the right hand side (rhs) of Eq. ([T]) represents the change in the 
average number of the vertices with degree k and fitness x due to process (i). The 
second term on the rhs accounts for the continuous introduction, with rate p, of 
new vertices with fitnesses drawn from the probability distribution f(x). The last 
term on the rhs represents the change in the average number of the vertices with 
degree k and fitness x due to process (ii). We also define 



k=0 



N(x,t)=J2N k (x,t), (2) 

and 

N(t) = I N(x,t)dx, (3) 



as the average number of the vertices of fitness x at time t and the average number 
of the vertices at time t, respectively. Summing Eq. ([!]) over k we obtain 



which yields 



dN(x,t) 

— gl — =Pf{x), (4) 



N(x,t)=pf(x)t + mf(x). (5) 



Integrating Eq. ([S]) we find as expected N(t) = pt + m, and therefore Eq. §5§ can 
be rewritten as 



N(x,t) = f(x)N(t). 



(6) 
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Now we can rewrite the integrals in the third term on the rhs of Eq. ([1} in terms 
of / and g 

^HZJ^ g{^x')N k {x\t)dx' 
EfcLo io°° Efelo Jo°° 9(x, x')N k (x', t)dx'N k (x, t)dx 

1 ^ g{x,x')f(x')dx' 

N(t) g(x,x>)f(x)f(x>)dxdx> 



N(ty 



(7) 



where 



A(x) = If 9(x,x')f(x')dx' (g) 
Jo°° Jo°° 9( x > x')f(x)f(x')dxdx' 



Furthermore, we assume that 



B(t)= f(x')J2 g(x,x')N k {x,t)dxdx', (9) 

^0 ,,_ n JO 



fc=0 

and the differential of which reads 
dB(t) 



P fix') g(x,x')f (x)dxdx' =C, (10) 



dt 

where C is a constant. Thus N k (x,t) grows linearly with t, and we introduce the 
function D k {x) such that 

N k (x,t) = D k (x)t. (11) 

The degree distribution of vertices can be estimated from D k {x) instead. Substitut- 
ing Eqs. ([7]) and (jlll) into Eq. ([T]) gives the recursive equation for D k (x) 



Ph,if(x) 



By defining 



(12) 



c 

we can rewrite Eq. ([121) as 



(x) = V - [°° f(x')g(x, x')dx' + ^-J?±A{x), (13) 



ff(x) . , pS k ,if(x) , . 

DiAx) = — — i/Jt-i X H , 14 

u 7 il(x) + 1 F(x) + 1' 17 



which can be solved recursively 



p.f{x)H k -\x) 
Dk{x) = [H(x) + l]* ■ (15) 



6, 2010 5:25 WSPC/INSTRUCTION FILE 0905IJMPC 



Mutual Selection in Network Evolution: the Role of the Intrinsic Fitness 5 

The result demonstrates that for every fixed x, the degree distribution of the gen- 
erated network should follows the right-shewed behavior. Moreover, the mutual 
selection rule presented here brings on the proportionality of the vertex degree to 
its fitness, which means that H(x) is an implicit function of k. Thus, given proper 
forms of the linking probability function g{xi, Xj), one can construct networks with 
power-law degree distributions. 
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Fig. 1. (color online) Degree distributions of vertices of the generated networks for different fitness 
distribution functions: uniform (a), exponential f(x) = e~ x (b), and power-law f(x) = x -3 (c). 
The linking probability function is g(xi,Xj) = XjXj. Each plot corresponds to one experiment of 
network generation with parameters N = 10 5 and m = 10. 
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To test above argument, we present computer simulations of the model, as shown 
in Fig.[T] We choose the simplest case g(xi, Xj) = XiXj and plot degree distributions 
for three kinds distribution functions of vertex fitnesses: uniform, exponential, and 
power-law. Even for this basic form of g, one can still notice the generalized power 
laws of the degree distribution in all cases, in agreement with analytical predictions. 

In summary, we have presented an simple model to justify the ubiquity of SF 
networks in nature, which results from the mutual selection rule based on a symmet- 
ric linking probability function g(xi,Xj) dependent on the affinity of the intrinsic 
fitnesses of the involved vertices, i and j. We found that it is always possible to 
find a proper form of g so that the generated network is SF in spite of the fitness 
distribution. In case that the values of vertex degrees are not available, we believe 
that the present model is relatively suitable. 

The authors acknowledge financial support from NSFC (No. 10805033) and 
STCSM (No. 08ZR1408000). This work is sponsored by the Innovation Founda- 
tion of Shanghai University. 
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